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ABSTRACT 

Aims. The chemistry, distribution and mass of the gas in the transitional disk around the 5 Myr old B9.5 V star HD 141569A are constrained. 
Methods. A quasi 2-dimensional (2D) chemistry code for photon dominated regions (PDR) is used to calculate the chemistry and gas 
temperatures in the disk. The calculations are performed for several gas distributions, PAH abundances and values of the total gas mass. The 
resulting CO 7 = 2-1 and 7 = 3-2 emission lines are computed with a 2D radiative transfer code and are compared to observations. 
Results. The CO abundance is very sensitive to the total disk mass because the disk is in a regime where self-shielding just sets in. The 
observed CO emission lines are best fit by a power-law gas distribution of 80 Me starting at 80 AU from the central star, indicating that there 
is some gas in the inner hole. Predictions are made for intensities of atomic fine-structure lines. [C i], which is the dominant form of carbon in 
large parts of the disk, is found to be a good alternative tracer of the gas mass. 
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1. Introduction 

Circumstellar disks are the natural by-product of star forma- 
tion, and once a protostar has formed t hey are crucial for the 
further accretion of matter onto the star dShu et alJl987l) . Over 
the course of its life a circumstellar disk loses its gas until a 
dusty debris disk remains in which the small dust grains are 
produced by collisions of planetesimals. The transformation 
from gas-rich accretio n disks to dusty debris disks occurs rel- 
atively quickly (Simon & Pratdll995l) . and as a result only a 
few objects in this transitional stage have been found. One ex- 
ample is the disk around HD 141569A, which is a pre-main 
sequence Herbig Ae/Be star with spectral type B9.5 V at a 
distance of 99 pc. It has a luminosity of 22 Lq (Mermet al 
2QQ4J. Its age is estimat ed at ~ 5 + 3 Myr (.Weinberger et al 



2000t iMerfn etaDl2004 and it has two M-type companions 
HD 141569B and C (Rossiter 'l943). The optically thin dust 
disk arou nd HD 141569A has been observed in infrared emis- 
sion (e. g.lSvIvester et alJl996l) and in scattered light a t near in- 
frared ( Augereau et al. 1999; Weinberger et al. 19991 and vis- 
ible wavelengths ( Mouillet et al. 2001: Clampin et al. 2003.) . 
Evidence for both large dust grains JBoccaletti et ^I,|2Q(^ ^|) and 
polycyclic ar omatic hv drocarbons (PAHs. see 'Svlveste r et alJ 
[l996: .Weinberger et alj.2004) has been found. Gas in the disk 



has also been detected i n the form of CO p ure rotational lines 
JZuckerm an et alJ Il995t iDent et alJ l2005l) and ro-vibrational 
lines (^Bri' ttahretalJ l2003V Thus, this disk provides a unique 
laboratory to study the effects of dust evolution on the gas. 

Observations of the dust in the HD 141569A disk show 
a complex morphology: there is a large inner hole in the 
dust distribution out to 150 AU from the central star and 
two dust rings at 185 AU and 325 AU dAug^ereau et alJI 19991: 



Weinberger et al. '1999; Mou illet et"aD l200l"; 'B occaletti et'all 
2003: Clampin et al 2003.) . The outer dust ring may be ex- 



Send ojfprint requests to: 
jonkheid@strw. leidenuniv.nl 



B. Jonkheid, 



e-mail: 



plained by tidal interactions of the disk with the M-type com- 
panions ( Augereau & Paoaloi zou 2004), while both rings may 
also have formed du e to hy drodvnamic drag forces on the dust 
grains JxTahr & Liii l2005h or interaction with a giant planet 
|Wvatt"2005). The origin of the large inner hole is still unclear 
but it could contain PAHs (Weinberger et al., in prepa ration). 

Th e sii igle dish CO ob servations by Zuckerr nan et alJ 
il995) and IPent etalJ (^OO?) in principle provide constraints 
on the mass of gas in the disk. Since this object is in a tran- 
sitional stage it should be losing its gas rapidly, and determin- 
ing the gas mass and gas/dust ratio may give insight into the 
processes responsible for the gas loss. Additionally, the total 
amount of gas in the disk has implications on the dynamics of 
the dust grains, and an estimate of the gas mass and its distribu- 
tion may help future dynamical simulations of the disk. In this 
paper a chemical model is presented which calculates molecu- 
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Fig. 1. Input parameters of the disk model. Panel a gives the surface densities for dust (solid line) and for three different distribu- 
tions of gas (dashed lines). Panels b, d and f give the corresponding gas number densities. Panel c gives the mean dust temperature 
(averaged over the grain size distribution), where the white line denotes the 100 K isotherm. Panel e gives the mass density of 
dust grains. 



lar abundances at each position in the disk as well as line inten- 
sities to co mpare with submillimeter observations, building on 
the workbvlKamp & van Zadelhofil(l200ll) . lKamp et al.l(l2003l) 
and 'jonkheid et al.'('2004V By varying the gas distribution and 
the total gas mass of the disk, a best fit to the CO observations 
is presented. Predictions for other gas tracers are made. 

2. Model 

2.1. Input 

The basis for the calculations presente d here is forme d by the 
1 H- 1 -dimensional model described in Llonkheid et all lEo04) . 
In this model the disk is divided into a series of vertical 1- 
dimensional structures, each of which was treated as a pho- 
ton dominated region (PDR) illuminated from above. The 
PDR code calculates the chemistry using the detailed ra- 
dia tive transfer of an d CO disso ciating Unes described 
by black & van Dishoeckl ( 11987.) and Ivan Dishoeck & Blackl 
lll988l). The elemental abundances used in the calculations are 



Table 1. Adopted gas-phase elemental abundances with respect 
to hydrogen. 



Element 


abundance 


C 


1.3 X 10"* 


O 


2.9 X 10"* 


Mg 


4.2 X 10-*^ 


S 


1.9 X 10-*^ 


Si 


8.0 X 10^*^ 


Fe 


4.3 X 10"' 


PAH 


1 - 2x 10-'"^ 



" The PAH abundance varies between models with different gas 
masses, since the total mass in PAHs is kept constant. 



displayed in Table ^ it is assumed that the gas in the disk is 
of interstellar origin instead of second generation gas such as 
may result from the evaporation of solid bodies. The thermal 
balance is solved taking into account thermal coupling between 
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Fig. 2. Depr oiected HST image of t he HD 141569 disk in scat- 
tered light bv lMouillet et all ilQQlh . 



gas and dust, heating through photoionization of polycychc hy- 
drocarbons (PAHs), photoelectric effect on large grains, and the 
formation and dissociation of H2 and cooling through CO ro- 
tational lines and the fine-structure lines of C, and O. The 
code was tested extensively against other codes in a benchmark 
project (Roellig et al., in preparation). 

Since the HD 141569A disk differs significantly in terms 
of mass and structur e from the T-Tauri disks examined in 
Ijonkheid et all tOOA a different input is used (see Fig ure [T). 
The dust distribution is taken from lMouillet et al](l200ll) bv de- 
projecting the disk assuming a scattering phase function with 
g - 0.2 in the V-band (see Figure |2ji. The dust content in the 
inner hole is uncertain: the upper limit to the surface density is 
2 X 10"* g cm"^. In the calculations a value of 2 x 10"^ g cm"^ 
is used at R < 150 AU to be on the conservative side. An az- 
imuthal average of the distribution was then taken to obtain the 
surface density in Figure [ij. 

Several trial profiles were used for the radial gas distri- 
bution: a power law similar to the minimal mass solar neb- 
ula where Sgas «: R~^'^, with an inner radius of either 15 AU 
(called distribution I, Figure [Q?) or 80 AU (distribution II, 
Figure [ijl), and a distribution similar to the dust distribution 
(with an effective inner radius of 150 AU, called distribution 
III, see Figure^O- For each assumption of the radial distribu- 
tion the total gas mass was varied by scaling the entire distri- 
bution by a single factor The densities of gas and dust were 
derived from their surface densities assuming a vertical distri- 
bution of n oc e ^ where the scale height h - 0.085 R and R 
and z are cylindrical coordinates. The densities shown in Figure 
[l]for the different radial distributions are for total gas masses 
of 80 M(g. The dust temperature is determined assuming the 
grains are in thermal balance and the disk is optically thin in 
continuum radiation at infrared through ultraviolet (UV) wave- 
lengths. 

The in terstellar spect rum used for the UV radi- 
ation by Ijonkheid et al ] 1I2OO4 is inappropriate for 
HD 141569 in the chemically important wavelength range of 
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Fig. 3. The flux at the inner edge of the gas disk in distribu- 
tion I and II (15 AU from the central star), for the radiation 
field of a B9.5 star (solid line) and of the interstellar radiation 
field (dotted line), scaled to fit the integrated flux between 912 
and 3600 A, corresponding to /uv = 10^ times the interstellar 
radiation field. The spectral regime where H2 and CO are dis- 
sociated and C ionized, and the regime where CH is dissociated 
are indicated with bars at the top. 

912A < A < 1100 A, s o the spectrum of an AO star (taken from 
'Hauschildt et al."l999', with T^s = 10000 K and log^ = 4, 
see Figure |3} was used here to calculate the photoionization 
and photodissociation rates. For the photoelectric heating 
rate the differences in the spectral shape of the radiation field 
were i mplemented using the correction factor by Spaans et alj 
( I1994I) (Equation (13) in that paper). The radiation field has a 
strength of /uv = lO^x the strength of the interstellar radiation 
field by Draine ( 1978) at 15 AU from the centi-al star; at 80 AU 
thi s factor is 3 x 10^ 

■Jonkh eid et ^ il2004h considered only more massive, opti- 
cally thick disks, so a number of changes had to be made in 
the code to simulate this transitional object. First, the vertical 
PDR structures used in the original code are inappropriate here 
because scattering will be far less important in an optically thin 
disk. Therefore the PDRs are now taken radially from the star, 
with the strength of the radiation field scaling with r^^ (where 
r is the distance to the star, not to be confused with R, the ra- 
dial component of the cylindrical coordinate system used in the 
plots) in addition to absorption effects. 

Second, the absorptions responsible for dissociation of H2 
and CO and the ionization of C may become optically thick 
if the column densities of these species toward the central star 
become sufficiently high. The interstellar radiation field then 
becomes the most important source of dissociation of these 
molecules. This is implemented in the model by taking the pho- 
todissociation rates for these molecules to be: 



R 



ph, tot 



= R 



ph, star 



+ R 



ph.ISM 



where /?ph,tot is the total photorate, /?ph,star is the photorate due 
to stellar light and i?ph,isM is the photorate due to the inter- 
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Fig. 4. Abundances (per H atom) for H2 (left column) C (middle column) and CO (right column) for models I (first row), II 
(second row). III (third row) and II without PAHs (fourth row). In all cases the gas mass is 80 Mg. 



stellar radiation field. The former rates are calculated using 
the detailed treatments of Black & van Dishoeck ( 1987) for H?, 
Ivan Dishoeck & Blackl (il988) for CO and van Disho eck ( 198^ 
for C. The latter rates are calculated assuming they have a value 
of 50% of the unshielded interstellar photodissociation rate at 
the disk's surface, with shielding calculated using the vertical 
column densities of H2, CO and C to ward the nearest surface as 
inputs f or the shielding function s bvlDraine & Bertoldil(ll996l) 
for H 7, Ivan Dishoeck & Blackl ( 1 19881) for CO and 'Werner 
ll 19701) for C. For all other processes the interstellar radia- 
tion field was ignored, because the stellar radiation fiel d has a 
strength of /uv > 1000 x the interstellar UV field of 'Draine" 
ifl978.) (disregarding absorption) even at 500 AU and domi- 
nates the continuum flux. 



2.2. Treatment of dust 

Based on th e sca ttering properties investigated by 
iBoccaletti etaP ( l2003l) and lAugereau & P apaloizou {200^ 
and the SED fitting bv lLi & Lunind ll2003l) the dust particles in 
the HD 141569A disk are found to be much larger than their 
counterparts in the interstellar medium: HD 141569 dust has 
grain sizes 1 fim < a < I cm (with a the equiva lent spherical 
radius), and a size distribution n(d) <x a^^^ jLi & Lunind 
2003), while the ISM grains ha ye sizes 5 nm < a < 0.25 fim 
and a distribution n(a) oc a^^^ (tMathis et al.ll 19771) . Since the 
Leiden FDR code assumes dust grains to have an interstellar 
size distribution, it has to be modified to incorporate the larger 
sized grains. Dust grains enter the model in four ways: as 
absorbers for UV radiation, as sites for H2 formation, and in 
the thermal balance as a source for the photoelectric heating 
and gas-dust collisions. These processes are, to lowest order, 
proportional to either the surface area or the geometrical 
cross-section of the grains, so they can be approximated by 
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scaling them with 

\ /HD 141569 ' \ /ISM 

where < a"' >- j"'""" a'"n(a)da/ n(a)da, and assuming 
the number of grains to remain constant in a given volume. If 
the mass in grains is kept constant rather than the number, all 
rates have to be mutiplied with a further factor of 

»%rain,lSM ^ (""^)lSM 
'«grain,HD 141569 (1 " {^^}hD 141569 

where P is the porosity of the grains in the HD 141569A 
disk. Following Li & Lunine (2003), P = 0.73 for ice-covered 
grains. Since most of the dust has a temperature lower than 100 
K, this value is used everywhere in the disk. Thus, to lowest or- 
der all processes should be scaled by a factor 1/1500 using the 
grain parameters derived for HD 141569 for a gas/ dust mass 
ratio of 100. 

This treatment of the dust properties leaves some inaccu- 
racies of order unity, the most important ones being in the ex- 
tinction of stellar radiation and the photoelectric heating rate. 
The mean extinction cross sections and scattering properties 
were calculated in dependently using Mie theory with the di- 
electric function of iLi & Lunind ll200 3) for porous icy aggre- 
gates (Equation (9) in that paper); the results were used to cor- 
rect the extinction. The photoelectric heating rate was calcu- 
lated using a photoelectric yield of Y = 0.05 electrons/photon 
to simulate larger grains (Watso nll972h . For a single grain size 
of 3 um the resulting heating r ates are similar to those used by 
iKamp & van ZadelhoffI (l200ll) : for the HD 141569 grains the 
appropriate size distribution was used. 

Even though the mean dust temperature shown in Figure 
[l]is high enough to prevent freeze-out of CO on the grains, 
there remains the possibility that the large grains, which have 
a temperature lower than the mean dust temperature, act as 
sites for CO freezing. Calculations show that the temperature 
of > 100//m sized grains falls below 30 K only in the outer 
regions of the disk, and that even millimeter sized grains have 
temperatures > 25 K throughout the disk. It is therefore un- 
likely that freeze-out of CO will occur at any significant rate, 
so it is ignored in the calculations. 

From infrared observations by ISvlvester et alJ (Il996}) it is 
known that polycyclic hydrocarbons (PAHs) are present in the 
disk. In our model PAHs are assumed to be well mixed with 
the gas, i.e. the abundance of PAH molecules with respect to 
hydrogen is independent of location in the disk, even in the in- 
ner hole. The t otal amount of PAH s is kept constant at wjpah = 
7.9 X 10"*' Me jLi & Lunin3l2003l) . regardless of the total gas 
mass. This means that the PAH abundance varies between mod- 
els with different mass; their abundance is xpah = 1.5 x 10"'" 
per H atom when Mgas - 80 M^. PAHs affect the disk's chem- 
istry in three ways (see Jonkheid et al. 2004): as sources of pho- 
toelectric heating, as absorbers of UV radiation and by direct 
chemical reactions (particularly in the formation of H2 and the 
charge balance). These processes are not only important in the 
inner disk where practically no large grains are present, but also 
in the outer disk due to the decreased efficiency of the large 
grains. 
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3. Results 

3.1. Chemistry 

The results for the chemistry for models with different radial 
gas distributions are shown in Figures IHandlH It can be seen 
that models I and II show similar behaviour, with molecular 
abundances that are high at the inner edge of the disk. Model III 
on the other hand produces H2 and CO only very far in the outer 
regions of the disk. In all models the distribution of atomic car- 
bon resembles the H2 distribution. 

It can be seen in Figure|5lthat in most models the carbon in 
the midplane is mostly atomic rather than ionized or molecular 
(in the form of CO). The reason for this is that even though 
the radiation responsible for the photoionization of C and the 
photodissociation of CO is largely shielded by the upper layers, 
there are still many photons with A > 1200 A left over (see 
Figure|3}- These photons can proceed to dissociate precursors 
to CO (most notably CH), thereby decreasing the CO formation 
rate and thus driving the carbon chemistry to a neutral atomic 
phase. 

In Figure |5] the CO abundances in the midplane are also 
shown for model II with total gas masses of 60, 80 and 100 M^. 
It can be seen that the CO abundance depends strongly on the 
gas mass, since it changes by an order of magnitude with only 
a factor of 2 change in mass. This strong dependence on mass 
is due to self-shielding; in this regime shielding of CO just sets 
in, and therefore every small increase in CO abundance in the 
upper layers enhances the production of CO in lower layers, 
producing a positive feedback loop. 

Throughout the disk H2 is primarily formed on PAHs, and 
the large dust grains are very inefficient absorbers of UV ra- 
diation. The chemical structure therefore shows no correlation 
with the dust distribution. If PAHs are removed from the chem- 
istry H2 is primarily formed on large grains, and its abundance 
resembles the dust distribution. Because of the decreased effi- 
ciency of H2 formation on the large grains the molecular abun- 
dances are lower than when PAHs are present. This is espe- 
cially true for CO, which now cannot build up sufficient column 
densities to become self-shielding. 

It can be seen from Figure |6l that the contribution of the 
interstellar UV field to the photodissociation of H2 and CO 
and the photoionization of C can be more than a factor of 
100 higher than the stellar contribution in the outer disk, and 
therefore cannot be ignored in chemical calculations. This is 
especially true for C ionization, since the high abundance of 
atomic carbon in the inner disk shields the outer disk from the 
stellar radiation. The abundance of H2 is high enough for this 
molecule to become completely self-shielding in the vertical 
direction as well as the radial direction. In the case of CO and 
C, the vertical self-shielding only just sets in. Thus the relative 
contribution of the interstellar UV field to H2 dissociation is 
not as large as it is for CO dissociation and C ionization. 

The overall chemical timescales, which are controlled by 
the photorates, are short, ranging from from < 100 yr in the 
inner disk and surface layers to a few xlO^ yr in the midplane 
of the outer disk. Thus chemical equilibrium is attained on dy- 
namical timescales. 
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Fig. 5. Midplane abundances (per H atom) for H2 (upper left), CO (upper right) and C (lower left) for 80 M© gas distributions I 
(dotted line), II (solid line) and III (dashed line), and for distribution II when PAHs are removed from the chemistry (dash-dotted 
line). The lower right panel gives the midplane abundance of CO for gas distribution II with total gas masses of 60 (dotted line), 
80 (solid line) and 100 M© (dashed line). 



3.2. Temperature 

The results for the gas temperature are shown in Figure For 
models I and II the surface layers of the disk are cool (~ 20 K), 
because the photoelectric heating by PAHs and large grains is 
very inefficient here. In the lower layers the abundances of C 
and H2 increase, and the gas is heated by the photoionization 
and photodissociation of these species, thus increasing the tem- 
perature with respect to the surface layers (~ 50 K). In model 
III the structure is slightly more complex: in the inner parts of 
the disk, the cooling is inefficient due to the low densities there. 
In the outer disk the heating by photodissociation of H2 is very 
effective due to a combination of high abundances and high 
dissociation rate. It can be seen from Figure that removal 
of the PAHs from the thermal balance has little effect on the 
temperature structure. The temperatures found in the disk's in- 
ner edge match the CO rotational t emperature of 1 90 K found 
from the CO vibrational bands by Britta in et alJ (|2003). The 
temperatures found her e and their distribution agree well with 
the results of lKamp&'van Zadelhoft ( 2001.) for the Vega disk. 



Initially no temperature solution could be found in the in- 
ner disk in some models. It was found that the cooling rates 
for C^ and O were overest imated when the formulation of 
iTielens & HoUenbachl ( fl985l) was used: 

^x(vij) = riiAij hvij/3e^c (jij) C(vij) 

where A v is the cooling rate by species x, n, is the population 
density of species x in level /, A is the spontaneous transition 
probability, ySgsc is the escape probability and C(v,y) is a correc- 
tion factor for the infrared background radiation from hot dust 
and the cosmic microwave background radiation. In the Leiden 
PDR code UV pumping of the fine structure levels of C^ and 
O is also included in calculating the populations, and the above 
formula becomes invalid for very high UV fields. In the inner 
disk, the UV pumping becomes the dominant process for pop- 
ulating these levels, and as a result the cooling rates were too 
large. This was solved by formulating the cooling rate as 

^xivij) = hvijinjkji - riikij) 
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Fig. 6. The total photodissociation rates of H2 (panel a) and CO (panel c) and the total photoionization rate of C (panel e) in 
model II, and the ratio of the contibutions of the interstellar and stellar radiation fields (panels b, d and f). Results for models 
with different gas masses follow similar qualitative trends. 



where nj and «, are the populations of the lower and upper lev- 
els, respectively, and kji and kij are the collisional excitation 
and de-excitation rates. A stable temperature solution could be 
reached using this method. In the absence of UV pumping this 
reduces to the previous formulation. 

4. Constraining the gas mass 

4.1. CO lines 

Figure |S] shows the predicted CO 7 = 2 — > 1 and 7 = 3 — > 2 
submillimeter lines for 5 different models: the best fitting gas 
masses for the radial gas distributions I, II and III (see Figure 



0, and two additional values of the gas mass for the pre- 
ferred radial distribution II. The emission line profiles were 
calculated with the 2-dim ensional radiative transfer code by 
■Hogerheiide & van der Tak (2000 ). For these calculations, a 
distance to the HD 141569 system of 99 pc wa s used, and the 
disk is assumed to have an inclination of 55° (Moui llet et alJ 
EoOD . The results are convolved with a beam of 1 1 arcsec for 
the 2 ^ 1 line (IRAM 30m at 230 GHz), and of 14 arcsec for 
the 3 ^ 2 line (JCMT 346 GHz). Of the former line only the 
peak intensity and FWHM are known; of the latter, the model 
results are compared directly to the observed spectrum. The 
computational grid used in the radiative transfer code consists 
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Fig. 7. The gas temperatures found for models I (upper left), II (upper right), III (lower left) and II without PAHs (lower right). 



of 100 cells in the radial direction and 11 cells in the vertical 
direction. 

It can be seen that the radial distribution of the gas mainly 
influences the widths of the lines investigated here. Distribution 
I gives line profiles that are too broad to fit the observations, 
while those of distribution III are too narrow. The intermediate 
distribution II fits the linewidths nicely. Using this radial distri- 
bution the emission lines were calculated for several values of 
the total gas mass. The strong dependence of the CO abundance 
on the gas mass shown in Figure0]is reflected in the emission 
lines, where the peak intensity changes with a factor of 4 for 
less than a factor of 2 in mass. It can be seen that the 80 Me 
model gives a good fit to the CO 7 = 2^1 and 7 = 3 — > 2 
emission lines. The model predicts peak line intensities of 0.82 
and 0.78 K for the 7 = 4 - 3 and 7 = 6 - 5 Unes in a 10" 
beam. The model cannot account for the sharp peak seen in the 
red part of the 7 = 3-2 line, however This peak is likely 
due to asymmetries in the disk structure, which can only be in- 
vestigated by a 3-dimensional disk model and a 3-D radiative 
transfer code. 

There are several uncertainties in the derived gas mass of 80 
M®: first, the elemental abundance of gas-phase carbon is un- 
certain by a factor of a few since the fraction of carbon locked 
in grains is not known. Also, changes in other elemental abun- 
dances, and uncertainties in the PAH distribution and chemical 
rate coefficients affect the CO chemistry and thus the inten- 
sity of the CO lines in relation to the gas mass. Second, the 
CO chemistry is very sensitive to the stellar radiation field. If 
the star has a higher flux in the 912-1100 Arange this would 
directly aff'ect the abundances of C and CO. Third, the assump- 
tion that H2 can form on PAHs is critical to the chemistry pre- 
sented here. If the formation of H2 is less eflicient than assumed 
here it would reduce the abundance of all molecules in the disk. 



both because H2 is is a progenitor molecule for many species 
and because of the shielding to UV photons it supplies. For 
comparison of r esults obtained with a comple tely independent 
chemical code (Kamp & van ZadelhofflEooi . an uncertamty 
in the gas mass of a factor of a few is estimated. 

The effect of disk inclination on the line profiles is limited. 
Spectra of the CO 3-2 line for distribution II with a gas mass of 
80 Me with inclinations of 50° and 60° are very similar to the 
spectrum shown in Figure|8]and show equally good agreement 
with the observed spectrum. Since a 5° change in inclination 
would indicate a morphological eccentricity of the disk of ~ 
0.45, the value of 55° by Mouiflet et al. (2001) is retained. 

The models presented here confi rm the presence of gas in 
the inner d isk previously noted by 'Zuckerman et al.l ( Il995l) . 
iBrittain et al . (2003) and Dent et al . (2005). Furthermore, the 
presence of vib rationally exci t ed CO at temperatures of 150- 
200 K found bv lfirittain etaP ('2003') is matched by model II; 
model I has too high temperatures in the inner disk, while 
model III only has CO in the outer disk where the radiation 
field is too low to excite the vibrational levels of CO and the 
temperatures are too low to match the observations. 

4.2. Other emission lines 

The predicted [C i] and [C 11] fine structure lines at 609, 370 and 
158 fj.m are shown in Figure|5] The intensities were calculated 
assuming radial distribution II, and assuming the beam is ex- 
actly the size of the disk (10 arcsec). Although the intensities 
of the [C i] lines change less dramatically with disk mass than 
the CO lines they are still a good -and probably more direct 
(see Figurel^- indicator of the total disk mass. The [Cn] line 
on the other hand depends only very weakly on the disk mass; 
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Fig. 8. CO 7 = 2-1 (panels a and c) and 7 = 3-2 (panels b and d) emission lines. Panels a and b give the best fit emission 
lines for the three different radial gas distributions shown in Figure^ The dotted line is for model I, the solid line model II, and 
the dashed line for model III. Panels c and d give the emission lines for model II with total gas masses of 60 (dotted line), 
80 Mm (solid Une) and 100 (dashed line). The red lines give the pe ak intensity and fu ll width half maximum observed by 
IZuckerman et al.l 1 1995b in the left panels, and the spectrum measured bv bent et al] ( l2005l) for the right panels. 



while there is more carbon in disks of higher mass, the fraction 
of ionized carbon is smaller. 

It should be noted that the C abundance and [C i] intensity 
are sensitive to the UV intensity at A < 1100 A. In particular, 
the C^/C ratio is increased and the [C i] intensity decreased if 
the UV intensity drops less steeply at the shortest wavelengths. 
In that case, the CO self-shielding can be maintained with an 
only slightly higher gas mass while the relative importance of C 
photoionization to CH photodissociation is decreased, leading 
to a lower C abundance. 

4.3. Gas/dust ratio 

The spatial dust distribution considered in this paper has a to- 
tal mass of ~ 2.2 for grains with sizes 1 fim < a < I cm. 
Together with the gas mass of 80 this would give an over- 
all gas/dust mass ratio of 36, which cannot be distinguished 
from the interstellar value of 100 given the uncertainties in the 



gas mass. This may indicate that gas and dust are dissipated on 
similar timescales. However, the upper limit of the grain size 
distribution is not well constrained, and the total dust mass de- 
pends strongly on this parameter Furthermore, it is uncertain 
what part (if any) of the dust population is debris (i.e. the prod- 
uct of collisions between meter-sized bodies) and what part is 
"first generation" aggregates. Thus the derived gas/dust ratio 
has only limited meaning. 

4.4. Second generation gas? 

There is a possibility that the gas in an evolved disk such as 
that around HD 141569A, is entirely hydrogen-poor, second 
generation gas, i.e. gas originating from evaporating solid bod- 
ies instead of the original interstellar cloud. The CO photodis- 
sociation timescale is short, ~ 10^ years (see § I3.1> . In order 
to reproduce the observed CO mass of 4 X 10"* Me, CO has 
to evaporate from comets (or similar solid bodies) at a rate of 
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Fig. 9. The C and fine structure lines for Mgas - 60 (dotted line), 80 M® (solid line) and 100 M® (dashed Une). 



4x lO^^Meyr-i. If one takes a mass of 10'** g (2 x 10"'" Me) 
and a CO mass fra ction of 10% as t ypical numbers for solar 
system comets (see lDelsemrn3ll985l) this translates into a re- 
quired evaporation rate of 2 x 10^ comets per year, which seems 
too high to be reahstic. 

5. Conclusions 

A model of the HD 141569A disk has been constructed which 
calculates the chemistry and gas temperature in a quasi 2- 
dimensional formalism and predicts the resulting emission 
lines. The main conclusions are as follows: 

- A gas distribution with a total mass of 80 Me starting at 
80 AU from the central star produces the best fit to the ob- 
served CO emission lines for this object. Further observa- 
tions of gas-phase lines will be needed to corroborate this 
result. The gas distribution can be tested through spatially 
resolved CO observations. The [Ci] fine-structure lines are 
good tracers of gas mass, but the [Cn] fine-structure line 
shows only a very weak variation with disk mass. 

- The CO chemistry of the outer disk is very sensitive to the 
total gas mass due to the coupled effects of high densities 
(which drives carbon to its molecular form) and shielding 
of the stellar and interstellar radiation fields. 

- The interstellar radiation field plays a significant role in the 
dissociation of H2 and CO, because the lines through which 
these molecules dissociate become optically thick toward 
the central star. 

- PAHs are found to be very important in the chemistry. 
Because the grains in this disk are large, the most efficient 
way to form H2 is via PAHs, and therefore they affect the 
formation of all molecules. 

- Carbon is found to be predominantly in neutral atomic form 
due a combination of efficient self-shielding of carbon pho- 



toionization and rapid photodissociation of chemical pre- 
cursors of CO by the stellar radiation field. This conclusion 
depends strongly on the shape of the stellar radiation field 
at^ < 1100 A. 

- The current model reproduces the C O 7 = 3 — » 2 line 
profile observed bv lPent et alJ ( l2005b quite well. The sharp 
peak at the red side of the spectrum is not reproduced due 
to the cylindrical symmetry inherent in the model. Also the 
model uses a azimuthally averaged dust distribution while 
the real distribution is noticeably sharper, which also ex- 
plains the shallower profile. The current model ignores the 
possibility of a clumpy gas distribution, which would also 
result in sharper profiles. 

The calculations of the gas mass of the HD 141569 disk pre- 
sented here indicate that transitional disks retain a fraction of 
their original gas content after a few Myr. The remaining gas 
is not enough to form giant planets, but mechanisms exist for 
planets to form on shorter timescales than the age of the disk 
(e.g. Poll ack etalJfT9 96). It is expected that the remaining gas 
has a strong influence on the dust dynamics in the disk. 
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